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With the rapid developments of the all-solid-state deep-ultraviolet (deep-UV) lasers, the good nonlinear 
optical (NLO) crystal applied in this spectral region is currently lacking. Here, we design two novel NLO 
carbonates KBeC0 3 F and RbAlC0 3 F 2 from the first-principles theory implemented in the molecular 
engineering expert system especially for NLO crystals. Both structurally stable crystals possess very large 
energy band gaps and optical anisotropy, so they would become the very promising deep-UV NLO crystals 
alternative to KBBF. Recent experimental results on MNC0 3 F (M = K, Rb, Cs; N = Ca, Sr, Ba) not only 
confirm our calculations, but also suggest that the synthesis of the KBeC0 3 F and RbAlC0 3 F 2 crystals is 
feasible. 

Deep -ultraviolet (deep-UV, X < 200 nm) nonlinear optical (NLO) crystals, as a key component of all-solid- 
state deep-UV lasers, have played important roles in many advanced scientific and technical areas 1 " 3 . 
Numerous attempts have been performed on the explorations of new deep-UV NLO materials with good 
performances 4 " 8 , and the dominant research field has been focused on the borates 9 " 11 . In particular, Potassium 
beryllium fluoroborate (KBe 2 B0 3 F 2 , KBBF) exhibits excellent NLO performance in the deep-UV region; till now 
KBBF is the sole NLO crystal that can practically generate the deep-UV lasers by direct SHG process 1 . However, 
the applications of KBBF are heavily hindered by the layering tendency in the single crystal growth processes 12 . 
Therefore, it is urgently demanded the emergence of new types deep-UV NLO materials. 

To be a good deep-UV NLO material, the following four criterions are commonly considered 1314 : (i) a wide UV 
transparency range down to the deep-UV region, corresponding to the large energy band gap and high damage 
threshold; (ii) a relatively large efficient second harmonic generation (SHG) coefficient (dy ^ 0.39 pm/V, d 36 of 
KH 2 P0 4 ); (iii) a relative large birefringence (An > 0.08) to achieve the phase matching condition in the deep-UV 
region; and (iv) good chemical stability and mechanical properties. According to the anionic group theory 1516 , the 
planar [B0 3 ] 3 " microscopic anionic groups have the dominant contribution to the macroscopic optical aniso- 
tropics in crystal, as in the KBBF case. Analogously, [C0 3 ] 2 ~ and [N0 3 ]~ anionic groups are expected to be the 
good NLO micro -structural units as well since they have the similar planar triangle structure with the 7i-con- 
jugated molecular orbitals which can produce the large birefringence and second-order susceptibility. The 
nitrates, however, are not considered as the NLO candidates for their hydrolysis. Thus, it is greatly desirable 
to explore the deep-UV NLO materials in carbonates. 

The rapid developments of scientific computational resources make it possible to predict the new advanced 
materials directly from the first-principles theory 17 " 19 , although search for the best candidate for a special property 
is still a major task 20 . In this work, aiming at the additional NLO candidates in the deep-UV region, we design two 
structurally stable carbonates KBeC0 3 F and RbAlC0 3 F 2 . The energy band gaps and the linear and nonlinear 
optical properties of both crystals are determined by ab initio calculations with high precision, which clearly 
demonstrate their promising applications as good optoelectronic functional materials in the deep-UV region. 

Results 

In order to search for the suitable deep-UV carbonates, we employed the anionic group theory combined with 
first-principles calculations which are implemented in the molecular engineering expert system especially for 
NLO crystals originated from our group 915 . After numerous efforts, we have theoretically discovered two carbo- 
nates, KBeC0 3 F and RbAlC0 3 F 2 , according to the following materials design considerations: (i) the candidates 
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Figure 1 | Unit cell of KBeC0 3 F (a) and RbAlC0 3 F 2 (b). 



are the fluoride carbonates since they can possess the larger energy 
band gaps compared with other carbonates 21 ; (ii) the A-site cations 
are the metal cations without unclosed d or/electrons because the d- 
d or /-/electronic transitions have negative influences to the energy 
band gap. Especially the alkaline cations and the lightweight metal 
cations in the Group IIA and IIIA in the Periodic Table are taking 
into account; (iii) all [C0 3 ] 2 " anionic groups are parallel (flat-lying) 
with respect to the overall structural layering. This structural feature 
has strong optical anisotropic responses to the incident light and 
exhibit large optical birefringence; and (iv) the orientation of all 
[C0 3 ] 2 " groups is parallel to each other. This arrangement is favor- 
able to the additive superposition of the microscopic second suscept- 
ibilities in the anionic groups, so produce large macroscopic SHG 
effects in crystal. 

The structures of KBeC0 3 F and RbAlC0 3 F 2 are plotted in Fig. 1 
(their crystallographic data see Table SI in the Supplementary 
Information). In KBeC0 3 F, the [Be0 2 F 2 ] tetrahedra are connected 
with the [C0 3 ] triangles by sharing their edges, forming the zero- 
dimensional [BeC0 3 F 2 ] microscopic structures isolated by the K + 
cations. In RbAlC0 3 F 2 , the [A10 3 F 2 ] trigonal bipyramids and 
[C0 3 ] triangles are alternately arranged in a trigonal pattern and 
connected via common O corners, generating a two-dimensional 
infinite [A1C0 3 F 2 ] layer parallel to the a-b plane. The Rb + cations 
are located between these layers to balance charge and also hold the 
layers together through the coordination with O and F anions. Both 
structures clearly exhibit the parallel flat-lying arrangement of the 
[C0 3 ] 2 " groups in crystals, so they are expected to have large SHG 
coefficients and birefringences. More importantly, the structural sta- 
bilities of both crystals are carefully verified by the first-principles 
methods 18 ' 19 . Figure 2 displays the phonon spectra and the total 
energy as a function of volume per atom in KBeC0 3 F and 
RbAlC0 3 F 2 . In the phonon spectra (Figure 2(a)) none of the ima- 
ginary phonon modes is observed, and in the curves of the total 
energy as a function of volume per atom (Figure 2(b)) there has a 
single minimum in the wide range of volume modification. Both 
evidences clearly demonstrate that the titled crystals are kinetically 
stable. 

The calculated linear and nonlinear optical properties for 
KBeC0 3 F and RbAlC0 3 F 2 are listed in Table 1, and the experimental 
results for KBBF are also shown as a comparison. It is clear that 
both KBeC0 3 F and RbAlC0 3 F 2 possess very large energy band 
gaps, 7.61 eV (—164 nm) and 8.2 leV (—152 nm), respectively. 
Meanwhile, the birefringences and SHG coefficients in both 



KBeC0 3 F and RbAlC0 3 F 2 crystals are larger than those in KBBF. 
The large birefringences (An = 0.1297 in KBeC0 3 F and An = 0.0998 
in RbAlC0 3 F 2 ) guarantee the achievement of the SHG phase-match- 
ing condition in both crystals down to their UV absorption edge, thus 
their shortest phase-matching wavelengths very approach (or even 
exceed) the corresponding wavelength in KBBF (—161 nm) 12 . It 
should be emphasized that RbAlC0 3 F 2 is optically uniaxial, which 
would be beneficial to its practical applications. The excellent optical 
properties in KBeC0 3 F and RbAlC0 3 F 2 demonstrate that both crys- 
tals are comparable to KBBF for the deep-UV SHG capabilities. 
Upon being obtained, these crystals would have wide applications 
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Figure 2 | Evidences of the structural stabilities in KBeC0 3 F and 
RbAlC0 3 F 2 from the first-principles studies, (a) The phonon spectra, 
(b) The total energy as a function of volume per atom. The energy 
minimum is set as the zero point. 
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Table 1 | Comparison of the linear and nonlinear optical properties in KBeCOaF, RDAICO3F2 and KBBF. E g is the energy band gap, n is the 
refractive index at the wavelength ~ 400 nm, An is the birefringence, Xp^ is the shortest phase-matching wavelength, and c/ #/ - is the SHG 
coefficient. For a uniaxial crystal, n x = n y = n Q and n z = n e 



Crystal 


^ g (eV) 


n x 


n y 




An 


X PM (nm) 




4/(pm/V) 


KBeC0 3 F 


7.61 


1 .5458 


1.5149 


1.4161 


0.1297 


164 


c/ 2 i 


= -0.94; c/ 22 = 0.70; 


















= 0.26; c/ 12 = 0.15 


RbAIC0 3 F 2 


8.21 


1 .4982 


1 .4982 


1.3984 


0.0998 


152 


c/,6 


= -c/ 22 = "0.69 


KBBF a 


8.26 


1.4915 


1.4915 


1 .4035 


0.0880 


161 


du 


= 0.47 



a Ref. 12. 
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Figure 3 | (a) (left) Partial density of states (PDOS) projected on the constituent atoms in MNCO3F, exemplified by the KSrC0 3 F crystal, (right) 
the PDOS zoomed-in to show the non-bonding states on oxygen atoms, (b) PDOS projected on the constituent atoms in KBeC0 3 F (left) and 
RbAlC0 3 F 2 (right). The blue and red curves represent the s and p orbitals, respectively, and the shadow areas represent the region occupied by 
non-bonding states on oxygen atoms. 



SCIENTIFIC REPORTS | 3 : 1366 | DOI: 1 0.1 038/srep01 366 



3 



as the optoelectronic functional materials in the deep-UV spectral 
region. 

Discussion 

To verify the reliability of our ab initio calculations on the optical 
properties in carbonates, we also applied our calculated methods on a 
family of recently synthesized NLO fluoride carbonates, MNC0 3 ¥ 
(M = K, Rb, Cs; N = Ca, Sr, Ba) 22 . The preliminary optical measures 
revealed that these compounds possess wide energy band gap (UV 
absorption edge < 200 nm) and qualitatively exhibit good capabil- 
ities for the UV harmonic generation 22 . However, due to the very 
small size of the samples the obtained experimental data are not 
enough to determine their application prospects in the deep-UV 
region. Our calculated energy band gaps and SHG coefficients match 
the available experimental results very well (see Table S2 and S3 in 
the Supplementary Information), which strongly prove the validity 
and high precision of the first-principles studies on the UV NLO 
carbonate crystals. Since the MiVC0 3 F series possess strong NLO 
effect and large birefringence, we predict that they are suitable to 
be good NLO crystals and excellent birefringent materials in the UV 
region. Nevertheless, these carbonates cannot be applied in the deep- 
UV region due to their relatively small energy band gaps (E g ~ 
6.3 eV). 

The detailed electronic structure analysis show that the small 
energy gaps in MNC0 3 ¥ are mainly due to the O 2p non-bonding 
states exclusively occupied at the valence band maximum, which 
have negligible small overlap with other electronic states and directly 
determine the energy band gap 23 (see Figure SI in the Supplementary 
Information). The energy spanning of these non-bonding states is as 
large as about 2.5 eV, and their complete elimination can signifi- 
cantly increase the energy band gap of crystals to more than 
8.3 eV (UV absorption edge < 150 nm). Compared MNCO3F with 
our proposed crystals, the substitution of Ca, Sr or Ba atoms with the 
lightweight metal cations in the Group IIA or IIIA, Be or Al, can 
effectively remove the non-bonding orbitals in carbonates, analogous 
to the borate cases such as beryllium in KBBF 24,25 and aluminum in 
BaAlB0 3 F (BABF) 26 " 28 . Indeed, Figure 3 clearly shows that the 
energy spanning of the non-bonding regions in KBeC0 3 F and 
RbAlC0 3 F 2 reduces about 1.1 eV and 2.3 eV, respectively, com- 
pared with that in MiVC0 3 F, exemplified by KSrC0 3 F. Therefore, 
the UV absorption edge of both KBeC0 3 F and Rb A1C0 3 F 2 are sig- 
nificantly blue-shifted. In addition, the similarity of the micro-struc- 
tural features in the MiVC0 3 F series and our proposed crystals 
suggests that it is feasible to synthesize KBeC0 3 F and RbAlC0 3 F 2 
in experiments. 

We surveyed throughout Inorganic Crystal Structure Database 
(ICSD) 29 , and found that only 14 fluoride carbonates whose energy 
band gaps may be large enough to transmit the UV radiation have 
been discovered and synthesized in experiments (see Table S4 in the 
Supplementary Information), in which seven are noncentrosym- 
metric and satisfy the SHG requirement. The very few report on 
UV NLO carbonates actually implies that the researches on this type 
of optoelectronic functional materials have been long-termly 
neglected and the broad developing space is available. We believe 
that our work would have great implications on the search and design 
of new NLO crystals in the deep-UV spectral region with the great 
help of advanced synthesis and single- crystal growth approaches. 

Methods 

First-principles computational methods and geometry optimizations. The 

first -principles calculations are performed by the plane-wave pseudopotential method 30 
implemented in the CASTEP package 31 based on the density functional theory (DFT) 32 . 
The ion-electron interactions are modeled by the optimized normal-conserving 
pseudopotentials for all elements 33 " 35 . The kinetic energy cutoffs of 900 eV and 
Monkhorst-Pack /c-point meshes 36 with the spanning of less than 0.03/A 3 in the 
Brillouin zone are chosen. The supercell volume and the atomic positions for KBeC0 3 F 
and RbAlC0 3 F 2 are fully optimized using the quasi-Newton method 37 . The 
convergence thresholds between optimization cycles for energy change, maximum 



force, maximum stress, and maximum displacement are set as 10~ 5 eV/atom, 0.03 eV/ 
A, 0.05 GPa, and 0.001 A, respectively. The above computational set ups are sufficiently 
accurate for present purposes 38 . 

Calculated methods for structural stability. To verify the structural stability, 
phonon spectra calculations on KBeC0 3 F and RbAlC0 3 F 2 are performed. The 
positive eigenvalues for all phonon modes is the most important evidence for the 
structural stability in crystal. The structural stability can also be demonstrated by the 
total energy as a function of volume per atom. If there has only a single minimum in 
the wide range of volume modification, the studied geometrical structure would be 
stable. The studied methods on structural stability were also adopted by Yao 18 and 
Sheng 19 . 

Calculated methods for optical properties. When determining the linear optical 
properties the imaginary part of the dielectric function e 2 is calculated 39 , and then its 
real part is determined by Kramers -Kronig transform, from which the refractive 
indices (and the birefringence) are obtained. Moreover, the second-order 
susceptibility % (2) , i.e., the SHG coefficient dy, is calculated by the following formula 
developed by our group 40 : 

x *Py = x *Pv ( VE ) + x *ISy ( VH ) + % ^y ( tw0 bands ) ? 

where % a ^(VE), ^ aft '(VH) and % a/ ^(two bands) denote the contributions from 
virtual-electron processes, virtual-hole processes and two-band processes, respect- 
ively. The formulae for calculating x^VE), ^ a/3y (VH) and ^ a/?y (two bands) are given 
as follows 40 : 

**(vh)= x [S Wy)imK,pJy„] (-jV + -^—) 

n z m 5 z ^J4ti 3 or vc 

Here, a, /? and y are Cartesian components, v and v' denote valence bands, and c and c' 
denote conduction bands. P(cifiy) denotes full permutation and explicitly shows the 
Kleinman symmetry of the SHG coefficients. The band energy difference and 
momentum matrix elements are denotes as hajy andp|, respectively, and they are all 
implicitly k dependent. 

It is well acknowledged that the DFT calculations with the exchange-correlation 
(XC) functional of generalized gradient approximation (GGA) 41 always 
underestimate the energy band gap of crystals. For calculating the optical coefficients, 
a scissors operator 42,43 is usually introduced to shift up all the conduction bands to 
agree with the measured band gap, which is not "purely" ab initio. Recently, we have 
examined several XC functionals beyond GGA on the determination of energy band 
gaps and optical properties in UV NLO crystals 44 . It was found that the hybrid 
functionals such PBE0 45 , P3LYP 46 and sX-LDA 47 can predict the UV absorption edge 
very well, but the calculated electronic band structures and optical properties are not 
well reproduced compared to the scissors-corrected GGA method 44 . Indeed our 
calculations show that in the MATC0 3 F crystals the GGA energy band gaps are in large 
discrepancy to the experimental results (absolute error excesses 2.4 eV), while the 
PBE0 values are in good consistency with the measurements (see Table S2 in the 
Supplementary Information). Therefore, in this work the difference between PBE0 
and GGA energy band gaps is set as the scissors operator, which is then used to 
determine the optical properties by GGA. The further tests have also revealed that the 
X (2) values is insensitive to the modification (e.g. ~ ±0.5 eV) of the scissors operator 
because of the very large band gap (>6 eV) in the studied carbonates. Therefore, this 
computational procedure ensures that our first-principles studies are self- consistent 
without adjusting any parameter from experiments. 
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